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Abstract. By their very nature, field-theoretical Hamiltonians are derived in momentum representation. To 
solve the corresponding integro-difi'erential equations is more difiicult than to solve the simpler diflFerential 
equations in configuration space ('Schrodinger equation'). For the latter many diflterent and very efi'ective 
methods have been developed in the past. But rather than to Fourier-transform to configuration space - 
which is not always easy - the equations are solved here directly in momentum space, by using Gaussian 
quadratures. Special attention is given to the case where the potential in configuration space is linear 
and where the corresponding momentum-space kernel has an almost intractable l/(k — k')*-singularity. 
Its regularization requires a certain technical effort, introducing suitable counter terms. The method is 
numerically reliable and fast, faster than other methods in the literature. It should be useful to and also 
applicable in other approaches, including phenomenological Schrodinger-type equations. 

PACS. ll.lO.Ef Lagrangian and Hamiltonian approach - 11.15.Tk Other non-perturbative techniques - 
12. 38. Aw General properties of QCD (dynamics, confinement, etc.) - 12.38.Lg Other nonperturbative 
calculations 



1 Introduction 

Dirac's front form of Hamiltonian dynamics seems to be 
a good candidate for addressing to the thus far unsolved 
bound-state problem of a (gauge) field theory, particu- 
larly of non-Abelian quantum chromo-dynamics (QCD), 
as reviewed recently in 1^. Advantages are its compara- 
tively simple vacuum structure and the simple boost prop- 
erties. In fact, one can formulate the bound-state problem 
frame-independently, and if one uses the light-cone gauge 
A'^ = 0, the vacuum is trivial. More recently [||, the front- 
form Hamiltonian for QCD has been reduced to a renor- 
malized, effective Hamiltonian which acts only in the Fock 
space of one quark (g) and one anti-quark (q). One then 
faces the numerical problem to solve these equations. 

In principle, one could proceed hkc in for QED. 
But we prefer to move on slowly, by first suppressing in 
Section all fine and hyperfine interactions with a simple 
trick. The resulting spin-less equation finds its analogue in 
the 'central potential' of a Schrodinger equation. In Sec- 
tion ^, the so truncated effective interaction is shown to 
interpolate smoothly ~ as a function of only one numerical 
parameter ~ between a pure Coulomb and a linear poten- 
tial. These potentials in configuration space are presented 



in detail, together with the known analytical solutions to 
the Coulomb and the linear potential. Their knowledge is 
advantageous to check the numerical procedures. 

In Section ^ the counter term technology is presented 
in detail, and applied in Sections and ^ to the Coulomb 
and the Yukawa, and in Section to the combined prob- 
lem, respectively. The respective eigenvalues and eigen- 
functions are presented in all detail. More and perhaps re- 
dundant numerical results are compiled in the Appendix. 
A summary and a discussion in Section ^ rounds-off the 
paper. 



2 The light-cone momentum-space approach 

The front-form Hamiltonian for QCD has been reduced in 
to an effective Hamiltonian, with a resulting integral 
equation in momentum space 



+ 
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The spectrum of the invariant-mass squared eigenvalues 
is the goal we want tot reach. The corresponding 
eigenfunctions \'ip) are the probability amphtudes for find- 
ing a quark with longitudinal momentum fraction x, trans- 
versal momentum Icl , and helicity Xq , and correspondingly 
for the anti-quark. The physical (renormalized) masses of 
the quarks are denoted by rriq. The 4- momentum transfer 
along the quark line — —{k — k'Y{k — k')^ is generally 
different from Q?. Their mean = (Q^ + (5|)/2 appears 
in Eq.(]l|). The (renormalized) 'running' coupling constant 
a{Q) is a function of the momentum transfer and will be 
given below. The spinor factor 

(A,,A,|^(Q)|A;,A;-) (2) 

= [U{kq,\)^^^u{k'q,\'q)\ [v{kq, X,)^ ^v{k'q, X'^)] 

represents the current-current coupling and describes all 
fine and hyperfine interactions. 

To reduce the complexity of the problem we intro- 
duce several simplifications and approximations, as fol- 
lows. First, we suppress all fine and hyperfine interactions 
by setting 



[U{kq, Xq)ru{k'^A'q)\ [v{kq. ^.hM^'q, ^q)] 



=o(3) 



^ ^mqmqdx^\'^Sx-x'^ ■ 

The helicity summations in Eq.([l]) therefore drop out, and 
with equal quark masses rfiq = rfiq = m the equation 
simplifies to 



MV(a;,k_i) 



x{\ — x) 



^(x,kj_) 



(4) 



,/^^x)x'{l-x') 3~Q^ 



The longitudinal momentum fraction (0 < a; < 1) and the 
transversal momentum (— oo < k±^x < oo) have different 
domains of integration. It is convenient to transform 
integration variables (x,kj_) {kz^^l_) by 



(5) 



All components of the 'vector' k = (fc^, k^) have then the 
same domain. The corresponding Jacobian is 



dx 2 dk-. 



x{l — x) m A2(a;,kj_) 
with the A-factor defined by 



m^A*(a;,k_L) 



4x(l - x) 



+ k? + kt 



(6) 



(7) 



The wave function -0 can always be substituted by an other 
function 0, 



^{x,k^)^^l^^^{x,k^) . 

y^X{l ~ X) 



(8) 



Replacing the invariant mass-square eigenvalue by the 
binding energy E', i.e. 



4™ 



4mi? 



(9) 



the original integral equation Eq.(h|) becomes finally 



E 



1 



d^k' 



4 a(Q) 
27r2 J A(k)A(k') 3 g2 



(k') 



Note that the only approximation is the suppression of 
the fine and hyperfine interaction by Eq.(|^). 

At this point it is completely irrelevant that the equa- 
tion holds actually for the front form. One simply does 
not recognize its origin. It could be an equation in the in- 
stant form, where only the (usual) three momenta k play 
the role of integration variables. It would however not be 
trivial to Fourier transform this equation to configuration 
space, the factor A is preventing us to do that in a stan- 
dard way. We therefore apply one further approximation. 
Somewhat doubtfully, we apply the non-relativistic ap- 
proximation under the integral, and set 



A(k) ~ A(k') = 1, Q2 - (k - k')^ 



(11) 



One now has a complete analogy with the simpler spin- 
less case of an average potential. One could argue in fa- 
vor of this approximation, that it is only consistent with 
Eq.(||). By replacing the QCD running coupling 4/3 a{Q) 
with the (physical) QED coupling a one obtains the usual 
Coulomb-Schrodinger equation in momentum space. 

Amazingly enough, just by redefining the wave func- 
tions, we have derived an equation whose kernel is man- 
ifestly rotationally invariant. This is surprising since ro- 
tations perpendicular to the z-axis are dynamic operators 
in the front form and very complicated 

What shall be used for the effective coupling constant 
a{Q)l In an explicit expression was given, but here we 
use the simplified expression of Brodsky et al. M 



a(Q) 



127r 



127r 



25 In 



— — 



251n' 



25 In 



and adopt their numerical values 



(12) 



K ~ 160 MeV , ~ 872 MeV . (13) 

Now the problem is set. 

Approximating the logarithm In (1 -I- Q'^ / ji^) hy Q"^ / ^"^ 
we parameterize the kernel of integral equation (nG) as 

4a(Q) _ 



2s' 



(14) 



3 Q2 (k-k')2(c2 + (k-k')2) ■ 

With the parameter values of Eq. ( p^ ) one has 

s = 1070 MeV , c = 1605 MeV . (15) 



Note that for 



0, the coupling constant in Eq.(14) 



exposes a l/(k — k')^ singularity. 
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3 The potential in configuration space 

Before solving numerically the equation 

1 



m 



m = J 



k') (16) 



with 



we discuss its structure by transformation to configuration 
space. Applying the Fourier transformations according to 



^(x) = 
t/(k,k') = 



J d3qe''i^J7(k, k') 



(17) 



1 



7r2 (k-k')2((k-k')2+c2) 

one obtains a Schrodinger equation 

V2 - 



E 



2mr 



^A(x) = F(x)V(x) , 



(18) 



with the reduced mass rrir = m/2. The simple structure 
is a consequence of the kernel J7(k, k') depending only 
on the difference k — k'. Since / d^qe''^'' 1/(Q^ + c^) — 
— 27r2 e~'^''/r, with r = |x|, one has 



^(x) = ^ [ 



1] 



(19) 



where fi = 23^ jt? is similar to the usual fine-structure 
constant. The potential is a superposition of a Coulomb 
and a Yukawa potential, as visualized in Figs. |l| and ||. 
The behavior 



^(|x|) 



^ + s^r for r — > , 
-2s 1 for r — > oo 



(20) 



can be observed in the figures: For sufficiently small dis- 
tances the potential is linear (up to an additional con- 
stant), and for large r it becomes a Coulomb potential. For 
either of these extremes the analytic solutions are known. 



The eigenvalue E used in Eq. (|16|) depends on the three 
parameters to, s, and c. It is easy to show that it depends 
on them in the dimensionless combination 



2 _c 

fi TO 



By introducing the dimensionless variables 

fc , Q2 2E 



V 



we obtain 



TOr/3 ' 



to2/32 ' 



(21) 



(22) 



q2 _^ jy2 



cf{p') (23) 
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Fig. 1. The potential V{r) is plotted versus r for different 
values of 77 = 2c//3m.- Note that for r ^ the potential goes 
to the finite value — 2s^/c, with a linear slope. 
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Fig. 2. Same as Fig. g], but with the potential normalized to 
1/(0) = 0. 



which is much simpler than Eq. (y^ , and we will therefore 
focus on this version. 

It has the two solvable limiting cases: (1.) For r\ — > 
cx), the equation essentially shows the q^^ singularity of 
a Coulomb problem. (2.) For ry — > 0, the equation has a 
singularity corresponding to a linear potential, whose 
eigenf unctions are the Airy functions. 

The eigenvalues for the Coulomb problem are thus in 
our units 

e„ = ~4 ' " = 1,2,... . (24) 
The eigenvalues for the linear potential are 

e„ = -2?7 en?7^ , n = l,2,... . (25) 
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where ^„ are the zeros of the Airy functions Ai{£^n) = 0. 

The knowledge of these two limiting cases is very useful 
for testing the numerical results. 



4 Numerical solutions by Gaussian 
quadratures 

Eq.(p3|) is an integral equation in the three variables p. 
Restricting to s-waves one can integrate out the angles, 
which leads to an equation in one variable p = |p|, i.e. 



[e-p'] cj^ip)^- Hdp'?- 
Jo P 



Infc^ (26) 



iP + P') 



J\2 



In 



{p — p') + V 



{p ^p'Y + 772 



It is convenient to convert this integral equation for 
the unknown function 4>{j)) to a matrix equation for the 
unknown numbers Ui = ^JuJi4)i with (f)i = <j){pi), 



i.e. the transformation changes the weights uji into uji = 
c{yi)uJi. Thus 



Pi 



-z In ■ 



1 



and 



l + Vi 



(34) 



where the (hi and y,; are the tabulated weights and abscis- 
sas for the interval [—1, 1]. Finally Eq.(p9|) can be solved 
by conventional matrix diagonalization methods. 

Eq.(27) poses a problem: The diagonal matrix ele- 



ments an diverge logarithmically. Problems like these can 
be solved by the Nystr0m method - by the technique of 
counter terms - as to be discussed next ]lO| . In princi- 
ple one adds and subtracts in Eq.(|23|) a diagonal term 
F{p)(j){p). The idea is that one of them is treated analyt- 
ically and the other by Gaussian quadratures, such that 
the singularity in an cancels. 

The construction of suitable counter terms must be 
done separately for the Coulomb and the Yukawa part. 
Therefore we discuss these two cases in Sections || and ^ 
individually. In Section ^ we return to the full problem. 



z,j = l,...7V, (27) 



namely by approximating the integral with the finite sum 
of Gaussian quadratures. 



N 



f{p)dp — > ^Wi/(pO 



(28) 



The non-diagonal matrix elements are given by 



1 



Pj 



-^UJiUJj — 



In 



{Pi-Pj? 

{pi+pjy 



(29) 



In 



{pi +Pjy + tf_ 



It is numerically convenient to map the infinite interval of 
Eq.(p6|) onto a finite one by transforming variables ]lOt , 
according to 



f{p)dp 



f{p{y))'^dy 
1 ay 



(30) 



The mapping function yiji) is arbitrary, but must satisfy 
the boundary conditions 



y{p = 0) 



-1 



We choose 



yip) = 2e" 



y{p = 00) = -fl 



1 



(31) 



(32) 



with an adjustable 'stretching' parameter z. With cijj) 
dp/ dy, the quadratures then become 



^ f{pi)uji = ^ f{yi)c{yi)C 



(33) 



5 The Coulomb problem 

The Coulomb problem in momentum space was treated by 
Wolz as a numerical exercise and will be repeated 

briefly. In analogy to Eq.(E6h we want to solve 



p2] 0(p) = 



1 



dp' — In 

P 



ip-p') 



/^2^ 



ip + p')^ 



ct^ip'). (35) 



Adding and subtracting the analytically integrable Coul- 
omb counter term Fc{p)(t>{p) 



J (p - r 



(l+p2)2 



(1 



one arrives at 



{e-p^ + {i+p')) m = 
{p'-pf 



(36) 
(37) 

(38) 



7o P 



(p'+py 



4>{p') 



2^2 



(1+P'2)2 



<P{P) 



For p = p' , the term in the square bracket vanishes, i.e. it 
is justified to restrict to j ^ i in the summation over j on 
the r.h.s. of Eq.(p7|). The diagonal matrix elements then 
become 



TT ^ ^ 



Pj 1 (P 
LOj — In 

71" Pi 



2\2 



PjY (i+pf) 

(p, +p,)2 (l+p2)2 



(39) 



The off-diagonal matrix-elements in Eq.(g9|) are not af- 
fected by this procedure 



I_ rjj-^Pl_ (Pi ^Pjf 

* ^ Pi (Pi +Pjy 



(40) 



see Eq.(|29|). In the sequel we ask ourselves whether this 
results can be improved by a stretching factor z, as given 
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Exact 


Calculated 


n 


-Ijn' 


N=16 


N=32 


N=64 


N=16 


N=32 


N=64 






z = 0.70 


z = 1.0 


1 


-1.0000 


-1.0000 


-1.0000 


-1.0000 


-1.0000 


-1.0000 


-1.0000 


2 


-0.2500 


-0.2545 


-0.2525 


-0.2516 


-0.2523 


-0.2510 


-0.2506 


3 


-0.1111 


-0.1160 


-0.1135 


-0.1126 


-0.1145 


-0.1122 


-0.1117 


4 


-0.0625 


-0.0683 


-0.0649 


-0.0639 


-0.0676 


-0.0637 


-0.0630 


5 


-0.0400 


-0.0471 


-0.0425 


-0.0414 


-0.0485 


-0.0415 


-0.0406 



Table 1. The eigenvalues of the Coulomb problem for two values of the stretching factor z and three matrix dimensions A'^ are 
compared with the exact values. 




z 



Fig. 3. The Bohr eigenvalues for n = 2, 3, 4 are plotted ver- 
sus the stretching factor z for different matrix dimensions A'^ 

( N=8, N=16, N=32 and o N=64).- The exact 

eigenvalues —Xjr? are shown as well.- Note the hatched area, 
in which the numerical results are particularly stable against 
the number z. 



Fig. 4. The numerical eigenvalues for the Coulomb prob- 
lem are plotted versus the number of integration points 
(8, 16, 32, 64, 128, 256) for z = 0.70. 



the rate of convergence and the precision of the Gaussian 
method for solving equations in momentum space. 



in Eq.(34). The advantage of the stretching factor is that 
a particular choice shifts the bulk of integration points to 
the region where the wave function is significantly different 
from zero. In Fig. |3| the numerical eigenvalues of Eq.(^7|) 
with the matrix elements of Eqs.(p9|) and ( |4C| ) are plotted 
versus z for different matrix dimensions A^.^s seen in the 
figure, the functions e„(z) are rapidly varying (almost fluc- 
tuating) for low dimensionality and become flatter with 
increasing N . For a z within the hatched area, however, 
the numerical results are rather stable as function of N . 
A value of = 32 (or 16) seems to satisfy all practical 
requirements. The lowest eigenvalue (ei = —1.0) is not 
shown, since the function (ei(z)) is completely flat. These 
observations are somewhat more quantified in Table 0. 

The final results for the low lying part of the spectrum 
are accumulated in Fig. ^. The spectrum is remarkably 
insensitive to the matrix dimensions iV. Also the numerical 
wave function as displayed in Fig. ^is highly accurate. For 
all practical needs N — 16 is sufficient, see also Table |l|. 

One concludes that the particular choice of the stretch- 
ing factor (z = 0.70) can be important for improving 



6 The Yukawa potential 

The Yukawa potential in momentum space obeys the in- 
tegral equation 



[e-p \ (p[p) ^ - dp—\n - — ■ — — — - 



(41) 

In analogy to the Coulomb problem, one adds and sub- 
tracts an analytically integrable counter term Fy{p)4>{p) 



Fy ip) 



1 



1 



(1 + pT 



TT" J " (p - P')^ + (1 + P'^)^ 
(l+p2)2(p2 + (^_l)2) 



(42) 



~ (p2 -f 772 - 1)2 + 4p2 

The integral equation ( |4l| ) is then rewritten as 
{l+p^f{p^ + ij]-l)^)\ 



e-p 



(p2 7p 



4p2 
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Fig. 5. The exact wave function (j)(pi)\s for the Coulomb 
problem is plotted versus p (the momentum in units of the 
Bohr momentum), and compared with the numerical values 
Ui^/Aiujji (filled diamonds). Parameter values are z = 0.70 and 
A'' = 8.- Note the excellent agreement. 





4>{p''. 



dp'— In 
P 



ip' - P? + 



{p' + pY + "ff 

(l+p2)2 



(43) 



{i+p'^y 



The matrix equation (27) has then the diagonal elements 
„ _ (p| + (r]-l)2)(l+pf)2 



{Pl 



if - 1)2 +4p2 

{PJ~P^?+r|^^ i^+Pi? 



(44) 



{pj +PiY + rf 



(1 



-P^Y 



while the off-diagonal elements are not modified 

{P3-P^Y + rj^ 



1 



— ^ujiujj In 



{Pj +PiY + ri^ 



(45) 



The integral over the domain [0, oo] is mapped on the in- 
terval [—1, 1] as given in Eqs.(po|) to (^). The limit ?7 ^ 
gives back the Coulomb problem, see Section ^. We have 
checked explicitly that the programs reproduce this case. 
A comparatively small value of 77 is 77 = 0.01. Indeed, the 
low lying part of the spectrum and the wave function is 
very similar to the corresponding result for the Coulomb 
problem in Figs. ^ and ||. In either case, the eigenvalues are 
practically insensitive to iV, particular for the stretching 
parameter z = 0.7, chosen to be the same as for the Coul- 
omb problem. One should note that the Yukawa problem 
has only a finite number of bound states. The value of 
z — 1.25 is obtained by the same optimalization proce- 
dure and shown explicitly in Fig. ^. Finally we summarize 
our best values for 77 = 0.01 and 77 = 1.0 in Table ||. 



n 


77 = 0.01 


r)= 1.0 




z = 0.7 


z = 1.25 


1 


-0.8141 


-0.0208 


2 


-0.1022 


-0.0003 


3 


-0.0084 




4 


-0.0019 




5 


-0.0019 





Table 2. The eigenvalues of all bounded s-states for the Yu- 
kawa problem for two values of r). The matrix dimension is 
either iV = 32. 



0.01 - 



0.00 - 



-0.01 - 



-0.02 



-0.04 



-0.05 
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Fig. 6. The two lowest eigenvalues of the Yukawa problem 
with r] = 1.0 are plotted versus the stretching factor z for 

different matrix dimensions: ( N=8, ^_ N=16, N=32 

and o N=64). 



7 The hadronic Coulomb potential 

Combining appropriately the considerations for the Coul- 
omb and the Yukawa problem by adding and subtracting 
the counter terms {Fc{p) + FY{p))(j){p) as given in Sec- 
tions 1^ and ^, yield immediately the improved diagonal 
elements 



a^^ =-(l+p?)+ 277 + 

{i+p1? 



{p} + {ri-lf)(l+pW 



TTPl 



E 



ipl 

UjPj In 



f 772 - 1)2 

{Pi-Pj? 



^pi 



(46) 



In 



{pi - 
iP, 



2\2 



{I+P)) 



Pi? + 



[pj +PiY + v'^ 



■p'Y 



which together with the unchanged off-diagonal elements 
Uij as given in Eq.(p9|) define the matrix equation (p7|). 
In the limit 77 ^ 00 the spectrum of Eq. (p7|) is expected 



to be very close to the Coulomb spectrum, 



-IV 



fact, which has been used to test the computer codes. In 
Fig. it is shown that the Coulomb limit is already well 
achieved for numerical values of 77 > 2.0. 



S.Bielefeld et.al: On technically solving an effective Hamiltonian for QCD 



7 



20.0 



16.0 



12.0 



8.0 



/ 


/ 


Airy regime / 




4/3 / 










/ -1 +2 Tl 




Coulomb regime 


/ f 




/ y 





2.0 



4.0 



6.0 



8.0 



10.0 



12.0 



Fig. 7. The lowest eigenvalue t(r]) = e + 2?) is plotted versus r\ 
(o). The upper solid line indicates the Airy-type solution and 
the lower solid line visualizes the Coulomb-type solution. 



4.0 




Fig. 8. A zoom of Fig. ^ 



In the other limit, 0, the spectrum of Eqs.(p6[) 

or (^) is expected to approach the spectrum for a Unear 
potential, i.e. e„ — — 277 -I- ^„7y*/'^. A typical Airy-solution, 



< 



however, is achieved only for very small values, i.e. 
0.02, as seen in Fig. ||. 

For all other values, the spectrum is somehow interme- 
diate between those two extreme cases, as quantitatively 
demonstrated in Table ^. It is remarkable how the curve 
of the calculated eigenvalues e(ry) = 2ry -I- ei in Fig. ^ in- 
terpolates between the two asymptotic curves e(l) ^ Jy^^"^ 
and e(l) ~ 2??- 1. 

Note that the stretching factor z should be optimalized 
for each value of rj. The resulting value is given in Table I 
as well. In parenthesis we note that the Coulomb limit for 
the present hadronic case is reachable for 77 ^ 00 while 
for the Yukawa case it was 77 ^ (see Section ||). 



More explicit numerical results are given in App. 



8 Summary and discussion 

Since the components of total four-momentum commute 
with each other, a field theoretic Hamiltonian is formu- 
lated quite naturally in momentum representation. In the 
instant form (usual quantization) the constituent's de- 
grees of freedom are the three space-like momenta (and 
their helicities and flavors), in the front form (or in light- 
cone quantization) they are the longitudinal momentum 
fraction x and the two transversal momenta ki . Effective 
Hamiltonians have the same property, they become inte- 
gral equations in momentum space. Usually, one Fourier- 
transforms such a momentum-space integral equation to 
a Schrodinger-type equation in configuration space and 
solves it by the familiar methods. Taking Fourier trans- 
forms is however not always easy, if not impossible, with- 
out additional assumptions. 

In the present work we therefore want to solve the inte- 
gral equations directly in momentum space. In particular, 
we look at an integral equation with an interaction kernel 
like ^ 

C/(k,k') 



(k-k')2(c2 + (k-k')2)7r2 



as derived in Section ||. We want to calculate the eigenval- 
ues and eigenfunctions for the full range of the interaction 
parameters s and c, and of the physical mass m of the con- 
stituent particles, a quark and an anti-quark with equal 
mass. As shown in Section ^ the solutions are a function 
of only one dimensionless parameter 



which intern can be interpreted as the ratio of two di- 
mensionless parameters, c/m and s/c. The limit 77 — > 00 
corresponds to a pure-Coulomb kernel [/ ~ (k — k')"^, 
the limit 7/ — > generates the highly singular interaction 
kernel U ~ (k — k')~^, as typical for a linear potential. 
The latter case was also the topic of Hersbachs work Q . 

As demonstrated in several examples in Section ^ the 
transitional region with 77 ~ 1.0 corresponds to a superpo- 
sition of a Coulomb and a Yukawa potential, which both 
are studied on their own merit in Sections [s] and |^, re- 
spectively. Therefore, for values of s and c fixed by their 
default value in Eq. ( p^ ) , the spectra for a mass larger than 
typically 180 GeV are more like those in a linear poten- 
tial, as opposed to more Coulomb-like spectra for masses 
smaller than typically 1.8 GeV. Masses in between have 
led to a mixed pattern. 

The technical problem of solving the interaction ker- 
nel in momentum space is approached by discretization 
via Gaussian quadratures, and diagonalization of the so 
generated Hamiltonian matrix. Much attention is paid to 
speed up convergence by a counter-term technology al- 
ready developed for the pure Coulomb case . It is 
summarized in Sections ^ and |^. In Sections ^ and |^ it 
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V 


z 


e = e + 2?) 


Airy regime 


Coulomb regime 










-1 + 277 


0.001 


0.0370 


0.00021905 


0.00023381 




0.01 


0.0714 


0.00454474 


0.00503728 




0.1 


0.2174 


0.08860885 


0.10852499 




1.0 


0.3861 


1.48676105 


2.3381 


1 


2.0 


0.4762 


3.31000318 




3 


4.0 


0.8000 


7.16112492 




7 


10.0 


0.8333 


19.04647165 




19 



Table 3. The eigenvalues t -\- 2r) of the integral equation ( pfij ) are given for increasing analytical values of the physical 
parameter 77, Column 2 gives the actual stretching parameter z. In the last two columns the corresponding eigenvalue for the 
Airy or Coulomb regime are quoted for purpose of comparison. 



n 


N=16 




N=32 


N=64 


N=16 


N=32 


N=64 






0.01,z = 0.0714 




0.1, 2 = 0.2174 


1 


0.0045 




0.0045 


0.0046 


0.0884 


0.0886 


0.0887 


2 


0.0073 




0.0074 


0.0075 


0.1312 


0.1316 


0.1317 


3 


0.0092 




0.0094 


0.0095 


0.1541 


0.1548 


0.1549 


4 


0.0106 




0.0110 


0.0110 


0.1674 


0.1685 


0.1686 




r; = 


1.0, z = 0.3861 




10.0, z = 0.833 


1 


1.4860 




1.4868 


1.4871 


19.0456 


19.0465 


19.0470 


2 


1.8183 




1.8193 


1.8196 


19.7537 


19.7551 


19.7555 


3 


1.9090 




1.9102 


1.9105 


19.8877 


19.8898 


19.8903 


4 


1.9447 




1.9464 


1.9468 


19.9339 


19.9372 


19.9379 



Table 4. The spectrum e„ + 2?7 for different values of r\ and A^, at the optimalized value of 2. 



is adapted to the Yukawa and the combined Yukawa plus 
Coulomb problem, in this work refered to as the hadronic 
Coulomb problem by obvious reasons. Special emphasis is 
put on a free formal parameter, called the stretching fac- 
tor 0, which can be adjusted for a considerably increased 
numerical precision and stability. This way, one can re- 
strict - on the average ~ on matrix diagonalization prob- 
lems with a matrix dimension as small as iV = 32. This 
part of the technical problem is applicable to many other 
physical problems. The restriction to equal masses of the 
constituents can be relaxed easily. 

We find it remarkable, that we solve a problem which 
on the technical level looks like a problem of usual (equal- 
usual-time) quantization despite the fact that the gener- 
ated solutions hold for the front form. The equation actu- 
ally being solved cannot be recognized as to have its roots 
in light-cone quantization. The only physical assumption 
entering the considerations in Section g is that the current- 
current term was replaced by the term of leading order in 
Eq.(|). 

It is, of course, still a long way to go for solving an 
effective QCD Hamiltonian on the technical level. Only 
a few steps have been taken in the present work. On the 
long run, we want to include properly the non-local factors 
A(k) in Eq. (p^ , to relax the assumption of Eq.(||), and to 
insert more general expressions for the effective coupling 
constant a((3). Work in this direction is under way. 
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A Some selected numerical results 

In this appendix we add a few typical results for the ha- 
dronic Coulomb potential as discussed in Section |^. The 
spectra and s-wave eigenf unctions for a very small and 
a large value of 77 are shown to illustrate the main differ- 
ences. The structure of the spectrum for r\ = 0.01 is similar 
to those of a linear potential. With increasing values of 77, 
the shape of the spectrum changes to that of a Coulomb 
potential particularly for the large value r\ — 10.0. More 
results are available from the authors on request. 
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Fig. 9. The eigenvalues e + 2?; for 77 = 0.01 (z = 
0.0714) are plotted versus the number of integration points 
N (8,16,32,64,128,256).- Note the almost equidistant struc- 
ture as in a linear potential. 



Fig. 11. The eigenvalues e -I- 2?; for 77 = 10.0 {z = 0.833) are 
plotted versus the same number of integration points as left.- 
This spectrum is of the Coulomb type, see Fig. 
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Fig. 10. The s-wave eigenfunction for 77 — 0.01 is plotted 
versus the momenta in Bohr units.- Note that the most of the 
integration points are in that region where the wave function 
is significantly different from zero. This is due to the choice of 
z = 0.0714. 



Fig. 12. The s-wave eigenfunction for r\ = 10.0 (2 = 0.833) is 
plotted versus the momentum values p/pBohr- This figure is 
comparable with Fig.H. 
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Fig. 14. S-wave eigenfunction for r? = 0.01 
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Fig. 19. Spectrum for ri = 1.0 



100 



200 



300 



12.0 



10.0 



Fig. 17. Spectrum for r? = 0.01 
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Fig. 20. S-wave eigenfunction for r] = 1.0 



Fig. 18. S-wave eigenfunction for r] = 0.1 
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Fig. 21. Eigenvalues for r\ = 0.01 versus stretching z 
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Fig. 23. Eigenvalues for r] = 1.0 versus stretching z 
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Fig. 24. Eigenvalues for r] = 10.0 versus stretching z 
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Fig. 22. Eigenvalues for r] = 0.1 versus stretching z 



